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Abstract 

Multiplicity distributions of neutral and charged particles arising from 
squeezed coherent states are investigated. Projections onto global isospin 
states are considered. We show how a small amount of squeezing can sig- 
nificantly change the multiplicity distributions. The formalism is proposed 
to describe the phenomenological properties of neutral and charged particles 
anomalously produced in hadronic and nuclear collisions at very high energies. 
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I. INTRODUCTION 



There has been a long-standing interest in the study of the multiphcity distributions 
of coherent states modified by various dynamical constraints The recent interest in 

the disoriented chiral condensate (DCC) as a possible explanation of the Centauro-type 
events has further extended the use of coherent states in the study of unusual charge 
distributions 0. It has also been suggested HTyjTTj] that a natural quantum description of 
the DCC can be given in terms of the squeezed states |T2|-|T^]. In this paper we investigate 
the neutral and charge distributions of squeezed coherent states constrained to form global 
isospin scalar and vector states that can possibly arise in nuclear collisions at very high 
energies. Our emphasis will be on the multiplicity distributions of the charged and neutral 
pions. It is therefore a significantly more detailed investigation than would be required, if 
the goal is only to determine the average fraction of neutral particles produced. 

If charged pions are always produced in pairs to maintain local charge neutrality of the 
source medium, it is very natural to consider squeezed states, which are created from the 
vacuum by operators that are exponentiations of the quadratic forms of the creation and 
annihilation operators. Although the squeezed states are constructed in quantum optics to 
reduce noise in certain applications [p!2|-p!5[1, their mathematical structure has features that 
are useful in quantizing classical fields. Their relevance to hadron production in high-energy 
collisions was recognized long ago @]. Recent interest in the subject has been stimulated 
by the study of chiral phase transition, where the long wavelength pion mode is found in 



the linear a model to grow rapidly after quenching ||l6ip!7| . It need not be a DCC, but 
the interest in developing a quantum description of DCC has led to the realization of the 
relevance of the squeezed states |1TI| . 

In our view that is rooted in the phenomenology of hadronic and nuclear collisions, the 
likelihood of the existence of a long wavelength pion mode or the creation of a DCC shielded 
from the normal vacuum by a hot shell in a heavy-ion collision is very remote, no matter how 
ideal or contrived the collision process may be. While we can think of various arguments 
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against the conditions favorable for their occurrence, there has not been any serious effort 
made to present a phenomenologically reasonable scenario in the form of realistic modeling 
that can demonstrate the feasibility of their appearance in actual collisions. Thus the subject 
is driven so far mainly by the attractiveness of the theoretical idea and the alluring possibility 
of solving concrete classical equations. Until such time as when experiments suggest the 
relevance of disoriented chiral condensates to the anomalous production of charge and neutral 
particles, we shall regard the subject as an interesting speculation that does not preempt 
other possible approaches. 

We begin with the observation of two facts. One is the empirical fact that anomalous 
events like the Centauro have been recorded in cosmic-ray experiments , though not in the 
laboratory. The other is the lack of a theoretical framework to describe the unusual charge 
distributions independent of the chiral dynamics. It seems to us useful to have a language 
to convey quantitatively the characteristics of the distributions without a rigid dynamical 
basis so that one may have what amounts to a representation in much the same way that 
the negative binomial distributions have been used to describe multiplicity distributions in 
general 0,|TB|JT^. While there exist earlier attempts to describe exotic neutral and charge 



distributions in terms of the coherent states it is better for our purpose to introduce 
some parameters that can be adjusted to change the neutral-to-charge ratio. To that end 
we find that the squeezed coherent states are well endowed with features suitable for our 
needs. Moreover, the consideration of global isotopic spin turns out to be crucial, as in 
|Jl[|, since unlike quantum optics the squeezed coherent states produced in hadronic and 
nuclear collisions must have a limited range of possible initial isotopic spin states. We 
should emphasize, however, that the dynamical origin of the squeezed coherent states in 
such collisions is not an issue that must be addressed first; we use the squeezed coherent 
state here only as a mathematical framework to discuss the production of neutral and charge 
quanta that are locally neutral in a certain limit (squeezed vacuum). 

We give in Sec. II a brief review of the squeezed coherent state. Application to par- 
ticle production is considered in Sec. Ill, and the results on multiphcity distributions are 



presented in Sec. IV. 



II. SQUEEZED STATES 



In this section we summarize the basic properties of squeezed states [p!4- [15| , |20|l that are 
to be used in the following section. 

We start with the coherent state la), created from the vacuum by the displacement 
operator D{a): 



\a) = D{a)\0) 



(2.1) 



where 



(2.2) 



and a are creation and annihilation operators, and a is a complex number. From ( p.2| 
we have 



D''{a) aD{a) = a + a 



(2.3) 



whereupon we obtain the fundamental property that defined in ( |2.1| ), is an eigenstate 
of a 



a\a) = a\a) 



(2.4) 



with eigenvalue a. 

The squeeze operator S{ri) is defined by 



S{ri) = exp {ri*a'^ — 7]a)'^)/2 



where rj = re^'^^. A unitary transformation on a by 3(7]) yields 



S^{ri) a S{ri) = jj,a — ua' 



(2.5) 



(2.6) 



where 
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/i = coshr, 1/ = e*^^sinhr . (2-7) 
Inverting ( p.6|) and a similar one on gives 

a{rf) = S{r]) a S^T]) = + v a) . (2.8) 
The unitary transformation of S{ri) on the displacement operator yields 

S\7^)D{a)S{r^) = D{a{7^)) , (2.9) 

where 

airf) = na + pa* . (2.10) 

The squeeze state |a,77) is defined by 

|a,r/) = D(«)S(r/)|0) . (2.11) 

Evidently, when = 0, i.e. no squeezing, |a,0) is identical to the coherent state \a). Even 

for 7^ 0, the expectation value of a is still a, since we have, using (|2.3| ) and (p^) 



(a,77|a|a,r7) = (0|5^(r/)(a + a)^(r/)|0) = a . (2.12) 

However, |a, rf) is not an eigenstate of a for r/ 7^ 0; instead, it is an eigenstate of 0(77), defined 
in dp), i.e., 

a(r/)|a,r/) = 5(r7)a5t(r7)Z}(a)5(r/)|0) 

= S{vi)a\a{vi)) = a{vi)\a,'q) . (2.13) 

Inverting ( p.8| ) we have 

a = ^ia{r]) — va\r]) (2-14) 

so that 

a)a = i^^i^ + |z/p) a\ini)a{7]) - ^iva^ivif - iju*a{r]f + |z/p , (2.15) 
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where 



{r]), a^{ri) = 1 has been used. It then follows from ( p. 13 ) and ( 2.1CI| ) that 



(^a,ri\a'a\a,ri^ = lal"^ + . (2.16) 

Thus the expectation value of the number operator for the squeezed state differs from that 
for the pure coherent state by the additive term |z/p. 

III. MULTIPLICITY DISTRIBUTIONS FOR SQUEEZED ISOSPIN STATES 

We now generalize the squeezed states described in the previous section to the isospin 
space in a manner similar to how the coherent states are generalized in Ref. 1. The annihi- 
lation operator a is now a vector in the isospin space, and the corresponding eigenstate \d) 
is created from the vacuum by the operator D{a). The squeezed state is then 

|a,r/) =D{a)S{Tim . (3.1) 

The complex vector a = an has a direction in the isospin space, specified by the polar angles 
6 and 0. A squeezed state with a specific total isospin / and 2;-component Iz can then be 
expressed in terms of |a, rf) in the usual way 

\a,r^)ij^ = I dnf,Ylj^{e,<P)\a{e,<P),Ti) . (3.2) 



What we need to do is to determine the projection of this state to the eigenstates of the 
number operators specifying the neutral (uq) and charged (ric) particle multiplicities. 
Then the multiplicity distribution of a squeezed isospin state is 
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PninA^^^v) = {no,nc\a,ri)jj^ . (3.3) 

It will be parameterized by the variables a and rj. More precisely, the results will depend 
only on the magnitudes |a| and r, as one might surmise from ( |2.16| ). 

To proceed, it is necessary to determine first {uq, nda, rj), which in turn can follow from 
a knowledge of the scalar product between a squeezed state and a coherent state. Following 
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Yuen [jT2[, we first note that the differential operator representation of any polynomial 
M {a, a)^ of a and in the coherent-state representation is 



Applying this to aij]) for M (a, a^), we have with the help of ( ^.13 ) 



{0\M{a,a^)\^p) = M[f| + ^^J*]{0\^) 



(3.4) 



(3 d 



^p\a(r))\d,r]J = a{ri) (^I3\d,ri 

This is a differential equation on ((3\d,ri) that has the solution 
[P\d,7]) = exp 



P\d,r] 



(3.5) 



1 




~2 


/3 



2 



1 

-/ 



(3.6) 



The last term that depends only on d{ri) and d*{ri) can be determined by working with 

d{ri) d 



d,r]\a{r])\p 



2 dd*{r]) 



d,r]\f3 



(3.7) 



The required form for / {d{ri), <y*{ri)) is therefore 



1 



f {d{r]),d*{ri)) = --\a{r])\ + —a {r]) . 



(3.8) 



2' 2/i 

— * 

Since P represents both the neutral and charged sectors, P cos6 and P sinO, in the isospin 
space, respectively, we have the usual representation of a coherent state in the \no,nc) 
basis, i.e., 

(no,n,|/?) = ^j=== {cosOr (sin^)"^ eH^I /2$(0) , (3.9) 

where $(0) = exp [i{n+ — n^)(f)\ is a phase factor that depends on the difference between + 
and — charged-particle multiphcities. 
Using (PI), and (p|) on 



(3.10) 



no,nc 
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and the expansion 



J2 ^T^Hn„{zo)HnSZc 



(3.11) 



where = t cos6, = t sin6 and similarly for zq and Zc, Hn{z) being the Hermite polynomial, 
we have 



{no, nc\a, rj) = A [a{vi)\ H„o 



Hr, 



(3.12) 



where 



(no+nc)/2 



exp 



(3.13) 



Note that An^^nc [(^{v)] is independent of the isospin direction, whereas the arguments of Hn 
are, since ao{ri) = a{ri)cos9 and adr]) = a{ri)sm9. For the isospin projected state (p.2|), we 
obtain correspondingly 



(no, n^\a, r])j r = A^^^n, (a, /x, u) 54^;J&) , 



where 



(3.14) 



b = a{r])/,j2^ 



(3.15) 



(3.16) 



When used in (|3.3|), this gives us the multiplicity distributions that we seek. 

We have not given any special attention to the normalization of (J3\a{r]fj as a solution 
to the differential equation in (|3.5| ). That freedom will be used in normalizing PngnS'^yV) 
at the end. 

There is more freedom in the squeezed states than is needed in our problem. From 
and (|2.1CI| ) we recall that 



a{rj) = (acoshr + Q;*sinhr e*^'^) 



n 



(3.17) 



Let us set a, /i, and u all real so that 



a(?7) = Q;(/i + z/)n , (3.18) 

and 

From (|3.3| ) and (|3.14|) we have 



b=-^^j==^ , /i = Vl + . (3.19) 



where is a normahzation factor. 
The average multiphcity (n) is 



{n) = {no) + {n,) = + , (3.21) 



as required by ( p.l6|) , which is invariant under rotation in the isospin space. We shall find 



it more convenient in the following to use the squeeze parameter s, defined by 

s = (u/af (3.22) 

so that 

(n) = 0^(1 + s) . (3.23) 

In all instances that we shall consider below, the value of s will be small, so (n) will not be 
effected very much by squeezing. Nevertheless, we shall see that for certain isospin states 
the no and ric distributions can be significantly infiuenced by the small amount of squeezing. 

In the next section we shall calculate Pn^nc o^ily for -^^ = and 1, although there exist no 
obstacles in carrying out the computation at higher values of /. For pp annihilation / < 1 is 
all that is necessary. For pp inelastic collisions the isospin of the pion state can in principle 
be as large as / = 2, but we shall restrict our consideration here to only / < 1. It is tacitly 
assumed that the analysis considered here is only appropriate for the particles produced in 
the central region. Even for nuclear collisions we hope that / < 1 is sufficient. 

For / = 0, 1, we can write 



(3.24) 



Thus the integration in (|3.15|) requires n+ — r7,_ = Iz- For the 6 integration let us define 

JnL(b) = dxx\l - x'y/'H^,ibx)H^SbVT^) (3.25) 

so that B^^^^^{b) is proportional to it. Let G{ti,t2) be the generating function, from which 
J^'r'l-n can be obtained: 



(3.26) 



tl=t2=0 



Eq. (p.25|) implies 



^(^1,^2) = E r dxx\l-X^y/^Hr,,{bx)Hr,^ (bVT^) 



no-f-ric 



no, Tic 

e 



1 

oo 



e " E (2b)"+- ^i^^^/ ^ + fc + \^ + ^ + 2 



(3.27) 



where B{u, v) is the Euler beta function. The factor e„+fc is 1 if n + = even, if n + = 
odd; it is due to the symmetry of the integrand under x ^ —x. It can be shown after some 
algebraic manipulations that 



dti ^ ^ 



-1) 



q—n 



ql 



r^.-n(O) . 



(3.28) 



t=o iq-n)\ 

Since the order of p of Hp{z) must be even when z = 0, q — n must be even, so (—1)''"" = 1. 
Using this in ([CTI) and {^^) yields 



no ric 



4'o'n.(&) = E E(2&)"+™^no-n(0)i/„.-„.(0) 
n=0 m=0 

/n + A; + l m + / + 2\ 

-D ? T: ^n+k ■ 



no Uric' 



(3.29) 



We finally have from ( p. 201) 



N' 



4 1 + 



-(7iO+"c)/2 2 



(3.30) 
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where A^' is another normahzation constant. From the exponents in ( |3.24|) and ( p.25|) , we 
have 

/ = : k = 0,1 = 0, n+ = n_ . 

J = l,/^ = 0: k = 1,1 = 0, n+ = n^ . (3.31) 
/ = 1, = 1 : k = 0, 1= 1, n+ = n^ + l . 

Clearly, Uc = n+ + n_ is even (odd), when n+ = n_ {n^ + 1). When k = 0, n must be even, 
so no must also be even; for k = 1, uq is odd. Thus the evenness or oddness of uq and 
can be summarized as follows: 

(a) 1 = 0: nQ = even, ric = even . 

(b) I = 1, Iz = : uq = odd, Uc = even . (3.32) 

(c) I = 1, Iz = 1 : tlq = even, ric = odd . 

IV. RESULTS ON MULTIPLICITY DISTRIBUTIONS 

In this section we present the results of our calculation of P^^^I^Ja, s). Our objective is 
to provide some insight into the dependences on squeezing and on the isospin states. To 
that end we shall fix the average multiplicity iji) at some reasonable value, which we take 
to be 25. Thus a and s are constrained through ( p. 23 ). To demonstrate the sensitivity 



on squeezing for certain isospin states, we shall consider only two small values of s, viz., 
s = 0.02 and 0.04, which is sufficient to reveal the qualitative features of the s dependence. 
Other higher values of s can readily be considered, but are not examined here to avoid visual 
confusion in the results to be presented. 

In calculating P^'^^^i^a, s) we note that s, as defined in (|3.22|) , is independent of the sign 



of u. In ( p.l9| ) we have taken fi and u to be real, but it does not preclude u to be negative, 
corresponding to = 7r/2 in ( p. 7]) . Since the result will depend on the sign of z/, we shall use 
the notation s±, to stand for sgn(z/) = ±, while specifying the value of s. In the following 
the distributions of the neutral and charged sectors are given separately, where for fixed (n) 
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(4.1) 



(4.2) 



For / = 0, the uq distribution (s±) is insensitive to s±, but the (s±) is highly 
sensitive to s±. This is shown in Figs. 1 and 2 for the cases z/ > and z/ < 0, respectively. 
Clearly, is peaked at low no, while the peak of shifts to low ric with increasing 
s_, although stationary in s+. The behavior of P^^^ for squeezed states is very similar 
to the coherent case considered in There, the asymptotic behavior is /~^^^, where 
/ = nQ/iriQ + ric), (obtained by using the Stirling formula in analytical estimates). Thus 
that law is not greatly influenced by slight squeezing. 

For / = 1, = 0, > 0, P^'^ (s+) is very nearly the same as P^^° (-5+), even though the 
former is only for odd and the latter for even rZc, as shown in Fig. 3. Their shapes depend 
sensitively on s+. For / = 1, = 0, z/ < 0, the distributions are very different; P^'^ (s_) 
is not sensitive to s_, but P^^°(s_) is. That is shown in Fig. 4. Upon comparing these 
to P^(,_„^, it is evident that the distributions depend crucially on the isospin, even though 
= in both cases. 

Finally, for / = 1, = 1, P^'^ {s±) behaves very similarly as P°^ ("S±), with peaking at 
low riQ and independence on s±, as is clear in Figs. 5 and 6. P^'^ (s±) is also similar to 

(-^i)' former is only for odd nc, while the latter is for even 71^. Thus the /"^''^ 

behavior is restored again. 

It is relevant to ask whether the rit = + ric ^ oo limit should be independent of 
squeezing. There are two cases to consider: s — and s finite. The s dependence can 
be seen by examining (3.29). For infinitesimal s, h is large, so the term with n = no and 
m = ric is dominant. Then the large nt behavior is independent of the Hermite polynomials. 
In the rif oo limit the s dependence of (3.29) cancels that in (3.20), resulting in no net 
dependence on s. But if s is set at a finite value initially, all terms in the sum in (3.29) can 
be important, since Hp{0) diverges with increasing p. The large Ut limit of can then 
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be very different and -P^q|^^ can depend in s, even if s is not large. In both cases, if Ut is not 
infinite, tfiere would in general be some dependence on s, the details of which can only be 
seen from the numerical results. It is for this reason that our results depend on squeezing. 

To summarize, we observe that in all isospin states the ric distributions are most sensitive 
to squeezing, not just in the value of s, but especially in the sign of z/. That is to be contrasted 
from the peaking of the neutral distribution P„q at Hq = for both the / = state and 
the I = 1, Iz = 1 state, revealing the preference for small numbers of neutral pions for 
either sign of u. On the other hand, for / = 1, = the neutral to charge ratio is about 
1 when u > 0, but is greater than 1 when z/ < 0. The latter therefore leans toward the 
characteristics of the anti-Centauro events. By varying J, J^, and s±, one can obtain a wide 
variety of possible shapes of P„^ that could be used to fit the observed charged particle 
distribution in high-energy collision for any given average total multiplicity. In that sense 
its usefulness is rather similar to that of the negative binomial distribution, which has a k 
parameter that adjusts the width. It should be noted that neither the negative binomial nor 
the pure squeezed- state distributions can reproduce properly the behavior of the cumulants 
observed in the experiments [^,^, whereas the question about squeezed isospin states is 



so far open. Unlike the negative binomial distribution, we now have a great deal more 
freedom to accommodate very anomalous production processes, where the neutral to charge 
ratio need not be 1:2. Indeed, P„o can be peaked at = 0, which is a feature that may 
render this description relevant to the Centauro events. By varying the squeeze parameter, 
one can not only adjust the neutral to charge ratio, but also change the detailed hq and 
dependences. 

We have not made any proposal in this paper to suggest that the dynamical origin for the 
anomalous production of particles is related to the squeezed states. We have only presented 
a formalism for the phenomenological description of unusual P„p and P„^. If in the future 
such a formalism is found empirically to be well suited for the data, then perhaps one may 
be motivated to look deeper into the question of whether squeezed states are produced in 
certain collisions at high energy. 
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Figure Captions 

Fig. 1 Multiplicity distributions of neutral and charged particles for 7 = and u > 0. 
Open symbols are for neutral particles and full symbols are for charged particles. 
Circles (open and full) are for s — 0.02; triangles (open and full) are for s = 0.04. 

Fig. 2 Same as in Fig. 1 but for / = 0, < 0. 

Fig. 3 Same as in Fig. 1 but for / = 1, = 0, z/ > 0. (a) and (b) are for both uq and 
Tie, but (c) is for ric only to show the dependence on s. 

Fig. 4 Same as in Fig. 3 but for I = 1, Iz — 0, u < 0. 

Fig. 5 Same as in Fig. 1 but for / = 1, = 1, t/ > 0. 

Fig. 6 Same as in Fig. 1 but for / = 1, = 1, i/ < 0. 
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